AIRE relies on Z-DNA to flag gene targets for thymic T cell tolerization.

AIRE is an unconventional transcription factor that enhances the expression of thousands of genes in medullary thymic epithelial cells and promotes clonal deletion or phenotypic diversion of self-reactive T cells1-4. The biological logic of AIRE's target specificity remains largely unclear as, in contrast to many transcription factors, it does not bind to a particular DNA sequence motif. Here we implemented two orthogonal approaches to investigate AIRE's cis-regulatory mechanisms: construction of a convolutional neural network and leveraging natural genetic variation through analysis of F1 hybrid mice5. Both approaches nominated Z-DNA and NFE2-MAF as putative positive influences on AIRE's target choices. Genome-wide mapping studies revealed that Z-DNA-forming and NFE2L2-binding motifs were positively associated with the inherent ability of a gene's promoter to generate DNA double-stranded breaks, and promoters showing strong double-stranded break generation were more likely to enter a poised state with accessible chromatin and already-assembled transcriptional machinery. Consequently, AIRE preferentially targets genes with poised promoters. We propose a model in which Z-DNA anchors the AIRE-mediated transcriptional program by enhancing double-stranded break generation and promoter poising. Beyond resolving a long-standing mechanistic conundrum, these findings suggest routes for manipulating T cell tolerance.

Flipons and small RNAs accentuate the asymmetries of pervasive transcription by the reset and sequence-specific microcoding of promoter conformation.

The role of alternate DNA conformations such as Z-DNA in the regulation of transcription is currently underappreciated. These structures are encoded by sequences called flipons, many of which are enriched in promoter and enhancer regions. Through a change in their conformation, flipons provide a tunable mechanism to mechanically reset promoters for the next round of transcription. They act as actuators that capture and release energy to ensure that the turnover of the proteins at promoters is optimized to cell state. Likewise, the single-stranded DNA formed as flipons cycle facilitates the docking of RNAs that are able to microcode promoter conformations and canalize the pervasive transcription commonly observed in metazoan genomes. The strand-specific nature of the interaction between RNA and DNA likely accounts for the known asymmetry of epigenetic marks present on the histone tetramers that pair to form nucleosomes. The role of these supercoil-dependent processes in promoter choice and transcriptional interference is reviewed. The evolutionary implications are examined: the resilience and canalization of flipon-dependent gene regulation is contrasted with the rapid adaptation enabled by the spread of flipon repeats throughout the genome. Overall, the current findings underscore the important role of flipons in modulating the readout of genetic information and how little we know about their biology.

Searching for New Z-DNA/Z-RNA Binding Proteins Based on Structural Similarity to Experimentally Validated Zα Domain.

Z-DNA and Z-RNA are functionally important left-handed structures of nucleic acids, which play a significant role in several molecular and biological processes including DNA replication, gene expression regulation and viral nucleic acid sensing. Most proteins that have been proven to interact with Z-DNA/Z-RNA contain the so-called Zα domain, which is structurally well conserved. To date, only eight proteins with Zα domain have been described within a few organisms (including human, mouse, Danio rerio, Trypanosoma brucei and some viruses). Therefore, this paper aimed to search for new Z-DNA/Z-RNA binding proteins in the complete PDB structures database and from the AlphaFold2 protein models. A structure-based similarity search found 14 proteins with highly similar Zα domain structure in experimentally-defined proteins and 185 proteins with a putative Zα domain using the AlphaFold2 models. Structure-based alignment and molecular docking confirmed high functional conservation of amino acids involved in Z-DNA/Z-RNA, suggesting that Z-DNA/Z-RNA recognition may play an important role in a variety of cellular processes.

Nonalternating purine pyrimidine sequences can form stable left-handed DNA duplex by strong topological constraint.

In vivo, left-handed DNA duplex (usually refers to Z-DNA) is mainly formed in the region of DNA with alternating purine pyrimidine (APP) sequence and plays significant biological roles. It is well known that d(CG)n sequence can form Z-DNA most easily under negative supercoil conditions, but its essence has not been well clarified. The study on sequence dependence of Z-DNA stability is very difficult without modification or inducers. Here, by the strong topological constraint caused by hybridization of two complementary short circular ssDNAs, left-handed duplex part was generated for various sequences, and their characteristics were investigated by using gel-shift after binding to specific proteins, CD and Tm analysis, and restriction enzyme cleavage. Under the strong topological constraint, non-APP sequences can also form left-handed DNA duplex as stable as that of APP sequences. As compared with non-APP sequences, the thermal stability difference for APP sequences between Z-form and B-form is smaller, which may be the reason that Z-DNA forms preferentially for APP ones. This result can help us to understand why nature selected APP sequences to regulate gene expression by transient Z-DNA formation, as well as why polymer with chirality can usually form both duplexes with left- or right-handed helix.

Z-DNA as a Tool for Nuclease-Free DNA Methyltransferase Assay.

Methylcytosines in mammalian genomes are the main epigenetic molecular codes that switch off the repertoire of genes in cell-type and cell-stage dependent manners. DNA methyltransferases (DMT) are dedicated to managing the status of cytosine methylation. DNA methylation is not only critical in normal development, but it is also implicated in cancers, degeneration, and senescence. Thus, the chemicals to control DMT have been suggested as anticancer drugs by reprogramming the gene expression profile in malignant cells. Here, we report a new optical technique to characterize the activity of DMT and the effect of inhibitors, utilizing the methylation-sensitive B-Z transition of DNA without bisulfite conversion, methylation-sensing proteins, and polymerase chain reaction amplification. With the high sensitivity of single-molecule FRET, this method detects the event of DNA methylation in a single DNA molecule and circumvents the need for amplification steps, permitting direct interpretation. This method also responds to hemi-methylated DNA. Dispensing with methylation-sensitive nucleases, this method preserves the molecular integrity and methylation state of target molecules. Sparing methylation-sensing nucleases and antibodies helps to avoid errors introduced by the antibody's incomplete specificity or variable activity of nucleases. With this new method, we demonstrated the inhibitory effect of several natural bio-active compounds on DMT. All taken together, our method offers quantitative assays for DMT and DMT-related anticancer drugs.

The Binding of Monoclonal and Polyclonal Anti-Z-DNA Antibodies to DNA of Various Species Origin.

DNA is a polymeric macromolecule that can display a variety of backbone conformations. While the classical B-DNA is a right-handed double helix, Z-DNA is a left-handed helix with a zig-zag orientation. The Z conformation depends upon the base sequence, base modification and supercoiling and is considered to be transient. To determine whether the presence of Z-DNA can be detected immunochemically, the binding of monoclonal and polyclonal anti-Z-DNA antibodies to a panel of natural DNA antigens was assessed by an ELISA using brominated poly(dG-dC) as a control for Z-DNA. As these studies showed, among natural DNA tested (Micrococcus luteus, calf thymus, Escherichiacoli, salmon sperm, lambda phage), micrococcal (MC) DNA showed the highest binding with both anti-Z-DNA preparations, and E. coli DNA showed binding with the monoclonal anti-DNA preparation. The specificity for Z-DNA conformation in MC DNA was demonstrated by an inhibition binding assay. An algorithm to identify propensity to form Z-DNA indicated that DNA from Mycobacterium tuberculosis could form Z-DNA, a prediction confirmed by immunoassay. Together, these findings indicate that anti-Z-DNA antibodies can serve as probes for the presence of Z-DNA in DNA of various species origin and that the content of Z-DNA varies significantly among DNA sources.

The Simple Biology of Flipons and Condensates Enhances the Evolution of Complexity.

The classical genetic code maps nucleotide triplets to amino acids. The associated sequence composition is complex, representing many elaborations during evolution of form and function. Other genomic elements code for the expression and processing of RNA transcripts. However, over 50% of the human genome consists of widely dispersed repetitive sequences. Among these are simple sequence repeats (SSRs), representing a class of flipons, that under physiological conditions, form alternative nucleic acid conformations such as Z-DNA, G4 quartets, I-motifs, and triplexes. Proteins that bind in a structure-specific manner enable the seeding of condensates with the potential to regulate a wide range of biological processes. SSRs also encode the low complexity peptide repeats to patch condensates together, increasing the number of combinations possible. In situations where SSRs are transcribed, SSR-specific, single-stranded binding proteins may further impact condensate formation. Jointly, flipons and patches speed evolution by enhancing the functionality of condensates. Here, the focus is on the selection of SSR flipons and peptide patches that solve for survival under a wide range of environmental contexts, generating complexity with simple parts.

Special Issue: A, B and Z: The Structure, Function and Genetics of Z-DNA and Z-RNA.

It is now difficult to believe that a biological function for the left-handed Z-DNA and Z-RNA conformations was once controversial. The papers in this Special Issue, "Z-DNA and Z-RNA: from Physical Structure to Biological Function", are based on presentations at the ABZ2021 meeting that was held virtually on 19 May 2021 and provide evidence for several biological functions of these structures. The first of its kind, this international conference gathered over 200 scientists from many disciplines to specifically address progress in research involving Z-DNA and Z-RNA. These high-energy left-handed conformers of B-DNA and A-RNA are associated with biological functions and disease outcomes, as evidenced from both mouse and human genetic studies. These alternative structures, referred to as "flipons", form under physiological conditions, regulate type I interferon responses and induce necroptosis during viral infection. They can also stimulate genetic instability, resulting in adaptive evolution and diseases such as cancer. The meeting featured cutting-edge science that was, for the most part, unpublished. We plan for the ABZ meeting to reconvene in 2022.

Crystal structure of Z-DNA in complex with the polyamine putrescine and potassium cations at ultra-high resolution.

The X-ray crystal structure of the d(CGCGCG)2/putrescine(2+)/K+ complex has been determined at 0.60 Šresolution. Stereochemical restraints were used only for the putrescinium dication, and 23 bonds and 18 angles of the Z-DNA nucleotides with dual conformation. The N atoms of the putrescine(2+) dication form three direct hydrogen bonds with the N7_G atoms of three different Z-DNA molecules, plus three water-mediated hydrogen bonds with cytosine, guanine and phosphate acceptors. A unique potassium cation was also unambiguously identified in the structure, albeit at a ∼0.5 occupation site shared with a water molecule, providing the first example of such a complex with Z-DNA. The K+ cation has coordination number of eight and an irregular coordination sphere, formed by four water molecules and four O atoms from two phosphate groups of the Z-DNA, including ligands present at fractional occupancy. The structural disorder of the Z-DNA duplex is manifested by the presence of alternate conformations along the DNA backbone. Comparison of the position and interactions of putrescine(2+) in the present structure with other ultra-high-resolution structures of Z-DNA in complexes with Mn2+ and Zn2+ ions shows that the dicationic putrescinium moiety can effectively substitute these metal ions for stabilization of Z-type DNA duplexes. Furthermore, this comparison also suggests that the spermine(4+) tetracation has a higher affinity for Z-DNA than K+.

Comparative review on left-handed Z-DNA.

Being polymorphic, deoxyribonucleic acid is worthy of raise a variety of structure like right-handed B to left-handed Z conformation. In left-handed contour of DNA consecutive nucleotides substitute between syn-arrangement and anti-arrangement, through the chain. 2D gel electrophoresis comprising d(PCpG)n of topo isomers of a plasmid inserts d(pCpG)n, in this 'n' ranges among 8 to 21, indicate the change of B-Z DNA. The high denseness of salt is required for conversion of B configuration d(CG)n toward Z configuration. The rate of B to Z transition is measured by "Cytosine Analogues" and "Fluorescence Spectroscopy". h-ZαADAR1 that a Z-DNA's binding domain, binds and stabilizes one part in Z configuration and therefore the remaining half in B deoxyribonucleic acid configuration. At halfway point, it creates B-Z junction. "Stacking" is the main reason for the B-Z DNA junction construction. Upregulation of ADAM-12, related with Z-DNA is said to a cause for cancer, arthritis, and hypertrophy. Z-DNA forming sequence (ZFS) conjointly generates massive - scale deletion in cells from mammals.

Systematic Approach to Find the Global Minimum of Relaxation Dispersion Data for Protein-Induced B-Z Transition of DNA.

Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion spectroscopy is commonly used for quantifying conformational changes of protein in µs-to-ms timescale transitions. To elucidate the dynamics and mechanism of protein binding, parameters implementing CPMG relaxation dispersion results must be appropriately determined. Building an analytical model for multi-state transitions is particularly complex. In this study, we developed a new global search algorithm that incorporates a random search approach combined with a field-dependent global parameterization method. The robust inter-dependence of the parameters carrying out the global search for individual residues (GSIR) or the global search for total residues (GSTR) provides information on the global minimum of the conformational transition process of the Zα domain of human ADAR1 (hZαADAR1)-DNA complex. The global search results indicated that a α-helical segment of hZαADAR1 provided the main contribution to the three-state conformational changes of a hZαADAR1-DNA complex with a slow B-Z exchange process. The two global exchange rate constants, kex and kZB, were found to be 844 and 9.8 s-1, respectively, in agreement with two regimes of residue-dependent chemical shift differences-the "dominant oscillatory regime" and "semi-oscillatory regime". We anticipate that our global search approach will lead to the development of quantification methods for conformational changes not only in Z-DNA binding protein (ZBP) binding interactions but also in various protein binding processes.

Water structure around a left-handed Z-DNA fragment analyzed by cryo neutron crystallography.

Even in high-quality X-ray crystal structures of oligonucleotides determined at a resolution of 1 Å or higher, the orientations of first-shell water molecules remain unclear. We used cryo neutron crystallography to gain insight into the H-bonding patterns of water molecules around the left-handed Z-DNA duplex [d(CGCGCG)]2. The neutron density visualized at 1.5 Å resolution for the first time allows us to pinpoint the orientations of most of the water molecules directly contacting the DNA and of many second-shell waters. In particular, H-bond acceptor and donor patterns for water participating in prominent hydration motifs inside the minor groove, on the convex surface or bridging nucleobase and phosphate oxygen atoms are finally revealed. Several water molecules display entirely unexpected orientations. For example, a water molecule located at H-bonding distance from O6 keto oxygen atoms of two adjacent guanines directs both its deuterium atoms away from the keto groups. Exocyclic amino groups of guanine (N2) and cytosine (N4) unexpectedly stabilize waters H-bonded to O2 keto oxygens from adjacent cytosines and O6 keto oxygens from adjacent guanines, respectively. Our structure offers the most detailed view to date of DNA solvation in the solid-state undistorted by metal ions or polyamines.

A widespread pathway for substitution of adenine by diaminopurine in phage genomes.

DNA modifications vary in form and function but generally do not alter Watson-Crick base pairing. Diaminopurine (Z) is an exception because it completely replaces adenine and forms three hydrogen bonds with thymine in cyanophage S-2L genomic DNA. However, the biosynthesis, prevalence, and importance of Z genomes remain unexplored. Here, we report a multienzyme system that supports Z-genome synthesis. We identified dozens of globally widespread phages harboring such enzymes, and we further verified the Z genome in one of these phages, Acinetobacter phage SH-Ab 15497, by using liquid chromatography with ultraviolet and mass spectrometry. The Z genome endows phages with evolutionary advantages for evading the attack of host restriction enzymes, and the characterization of its biosynthetic pathway enables Z-DNA production on a large scale for a diverse range of applications.

Sequence-dependent cost for Z-form shapes the torsion-driven B-Z transition via close interplay of Z-DNA and DNA bubble.

Despite recent genome-wide investigations of functional DNA elements, the mechanistic details about their actions remain elusive. One intriguing possibility is that DNA sequences with special patterns play biological roles, adopting non-B-DNA conformations. Here we investigated dynamics of thymine-guanine (TG) repeats, microsatellite sequences and recurrently found in promoters, as well as cytosine-guanine (CG) repeats, best-known Z-DNA forming sequence, in the aspect of Z-DNA formation. We measured the energy barriers of the B-Z transition with those repeats and discovered the sequence-dependent penalty for Z-DNA generates distinctive thermodynamic and kinetic features in the torque-induced transition. Due to the higher torsional stress required for Z-form in TG repeats, a bubble could be induced more easily, suppressing Z-DNA induction, but facilitate the B-Z interconversion kinetically at the transition midpoint. Thus, the Z-form by TG repeats has advantages as a torsion buffer and bubble selector while the Z-form by CG repeats likely behaves as torsion absorber. Our statistical physics model supports quantitatively the populations of Z-DNA and reveals the pivotal roles of bubbles in state dynamics. All taken together, a quantitative picture for the transition was deduced within the close interplay among bubbles, plectonemes and Z-DNA.

Remodeling Chromatin Induces Z-DNA Conformation Detected through Fourier Transform Infrared Spectroscopy.

The SWI/SNF complex is a highly conserved chromatin remodeling complex and can hydrolyze ATP by its catalytic subunit BRG1 or BRM to reconstruct the chromatin. To investigate whether this ATP-dependent chromatin remodeling could affect the DNA conformation, we therefore regulated (knocked down or overexpressed) BRG1/BRM in the cells and applied Fourier transform infrared (FTIR) spectroscopy to probe DNA conformational changes. As a result, we found that BRG1/BRM was indeed associated with the DNA conformational changes, in which knockdown of BRG1/BRM reduced Z-DNA conformation, while overexpression of BRG1/BRM enhanced Z-DNA conformation. This Z-DNA conformational transformation was also verified using the Z-DNA-binding proteins. Therefore, this work has provided a direct analytical tool to probe Z-DNA transformation upon ATP-dependent chromatin remodeling.

Protein-induced B-Z transition of DNA duplex containing a 2'-OMe guanosine.

Structural transformation of the canonical right-handed helix, B-DNA, to the non-canonical left-handed helix, Z-DNA, can be induced by the Zα domain of the human RNA editing enzyme ADAR1 (hZαADAR1). To characterize the site-specific preferences of binding and structural changes in DNA containing the 2'-O-methyl guanosine derivative (mG), titration of the imino proton spectra and chemical shift perturbations were performed on hZαADAR1 upon binding to Z-DNA. The structural transition between B-Z conformation as the changing ratio between DNA and protein showed a binding affinity of the modified DNA onto the Z-DNA binding protein similar to wild-type DNA or RNA. The chemical shift perturbation results showed that the overall structure and environment of the modified DNA revealed DNA-like properties rather than RNA-like characteristics. Moreover, we found evidence for two distinct regimes, "Z-DNA Sensing" and "Modification Sensing", based on the site-specific chemical shift perturbation between the DNA (or RNA) binding complex and the modified DNA-hZαADAR1 complex. Thus, we propose that modification of the sugar backbone of DNA with 2'-O-methyl guanosine promotes the changes in the surrounding α3 helical structural segment as well as the non-perturbed feature of the ß-hairpin region.

Short but Weak: The Z-DNA Lone-Pair⋠⋠⋠π Conundrum Challenges Standard Carbon Van der Waals Radii.

Current interest in lone-pair⋅⋅⋅π (lp⋅⋅⋅π) interactions is gaining momentum in biochemistry and (supramolecular) chemistry. However, the physicochemical origin of the exceptionally short (ca. 2.8 Å) oxygen-to-nucleobase plane distances observed in prototypical Z-DNA CpG steps remains unclear. High-level quantum mechanical calculations, including SAPT2+3 interaction energy decompositions, demonstrate that lp⋅⋅⋅π contacts do not result from n→π* orbital overlaps but from weak dispersion and electrostatic interactions combined with stereochemical effects imposed by the locally strained structural context. They also suggest that the carbon van der Waals (vdW) radii, originally derived for sp3 carbons, should not be used for smaller sp2 carbons attached to electron-withdrawing groups. Using a more adapted carbon vdW radius results in these lp⋅⋅⋅π contacts being no longer of the sub-vdW type. These findings challenge the whole lp⋅⋅⋅π concept that refers to elusive orbital interactions that fail to explain short interatomic contact distances.

Dynamics Studies of DNA with Non-canonical Structure Using NMR Spectroscopy.

The non-canonical structures of nucleic acids are essential for their diverse functions during various biological processes. These non-canonical structures can undergo conformational exchange among multiple structural states. Data on their dynamics can illustrate conformational transitions that play important roles in folding, stability, and biological function. Here, we discuss several examples of the non-canonical structures of DNA focusing on their dynamic characterization by NMR spectroscopy: (1) G-quadruplex structures and their complexes with target proteins; (2) i-motif structures and their complexes with proteins; (3) triplex structures; (4) left-handed Z-DNAs and their complexes with various Z-DNA binding proteins. This review provides insight into how the dynamic features of non-canonical DNA structures contribute to essential biological processes.

Emerging roles for R-loop structures in the management of topological stress.

R-loop structures are a prevalent class of alternative non-B DNA structures that form during transcription upon invasion of the DNA template by the nascent RNA. R-loops form universally in the genomes of organisms ranging from bacteriophages, bacteria, and yeasts to plants and animals, including mammals. A growing body of work has linked these structures to both physiological and pathological processes, in particular to genome instability. The rising interest in R-loops is placing new emphasis on understanding the fundamental physicochemical forces driving their formation and stability. Pioneering work in Escherichia coli revealed that DNA topology, in particular negative DNA superhelicity, plays a key role in driving R-loops. A clear role for DNA sequence was later uncovered. Here, we review and synthesize available evidence on the roles of DNA sequence and DNA topology in controlling R-loop formation and stability. Factoring in recent developments in R-loop modeling and single-molecule profiling, we propose a coherent model accounting for the interplay between DNA sequence and DNA topology in driving R-loop structure formation. This model reveals R-loops in a new light as powerful and reversible topological stress relievers, an insight that significantly expands the repertoire of R-loops' potential biological roles under both normal and aberrant conditions.

A Genetic Instruction Code Based on DNA Conformation.

Flipons are sequences capable of forming either right- or left-handed DNA under physiological conditions, forming a class of dissipative structures that trade metabolic energy for information by cycling DNA between different chromatin states. Flipons enhance the diversity of transcriptomes, increasing entropy while enabling the evolution of features both new and unexpected.